A QM/MM approach for the study of monolayer-protected gold clusters

We report the development and implementation of hybrid methods that combine quantum mechanics (QM) with molecular mechanics (MM) to theoretically characterize thiolated gold clusters. We use, as training systems, structures such as Au25(SCH2-R)18 and Au38(SCH2-R)24, which can be readily compared with recent crystallographic data. We envision that such an approach will lead to an accurate description of key structural and electronic signatures at a fraction of the cost of a full quantum chemical treatment. As an example, we demonstrate that calculations of the 1H and 13C NMR shielding constants with our proposed QM/MM model maintain the qualitative features of a full DFT calculation, with an order-of-magnitude increase in computational efficiency.Comment: Journal of Materials Science, 201

Kinetics of diffusion-controlled enzymatic reactions with charged substrates

The Debye-Hückel limiting law (DHL) has often been used to estimate rate constants of diffusion-controlled reactions under different ionic strengths. Two main approximations are adopted in DHL: one is that the solution of the linearized Poisson-Boltzmann equation for a spherical cavity is used to estimate the excess electrostatic free energy of a solution; the other is that details of electrostatic interactions of the solutes are neglected. This makes DHL applicable only at low ionic strengths and dilute solutions (very low substrate/solute concentrations). We show in this work that through numerical solution of the Poisson-Nernst-Planck equations, diffusion-reaction processes can be studied at a variety of conditions including realistically concentrated solutions, high ionic strength, and certainly with non-equilibrium charge distributions. Reaction rate coefficients for the acetylcholine-acetylcholinesterase system are predicted to strongly depend on both ionic strength and substrate concentration. In particular, they increase considerably with increase of substrate concentrations at a fixed ionic strength, which is open to experimental testing. This phenomenon is also verified on a simple model, and is expected to be general for electrostatically attracting enzyme-substrate systems

An exact expression to calculate the derivatives of position-dependent observables in molecular simulations with flexible constraints

In this work, we introduce an algorithm to compute the derivatives of physical observables along the constrained subspace when flexible constraints are imposed on the system (i.e., constraints in which the hard coordinates are fixed to configuration-dependent values). The presented scheme is exact, it does not contain any tunable parameter, and it only requires the calculation and inversion of a sub-block of the Hessian matrix of second derivatives of the function through which the constraints are defined. We also present a practical application to the case in which the sought observables are the Euclidean coordinates of complex molecular systems, and the function whose minimization defines the constraints is the potential energy. Finally, and in order to validate the method, which, as far as we are aware, is the first of its kind in the literature, we compare it to the natural and straightforward finite-differences approach in three molecules of biological relevance: methanol, N-methyl-acetamide and a tri-glycine peptideComment: 13 pages, 8 figures, published versio

The regulatory subunit of PKA-I remains partially structured and undergoes β-aggregation upon thermal denaturation

Background: The regulatory subunit (R) of cAMP-dependent protein kinase (PKA) is a modular flexible protein that responds with large conformational changes to the binding of the effector cAMP. Considering its highly dynamic nature, the protein is rather stable. We studied the thermal denaturation of full-length RIα and a truncated RIα(92-381) that contains the tandem cyclic nucleotide binding (CNB) domains A and B. Methodology/Principal Findings: As revealed by circular dichroism (CD) and differential scanning calorimetry, both RIα proteins contain significant residual structure in the heat-denatured state. As evidenced by CD, the predominantly α-helical spectrum at 25°C with double negative peaks at 209 and 222 nm changes to a spectrum with a single negative peak at 212-216 nm, characteristic of β-structure. A similar α→β transition occurs at higher temperature in the presence of cAMP. Thioflavin T fluorescence and atomic force microscopy studies support the notion that the structural transition is associated with cross-β-intermolecular aggregation and formation of non-fibrillar oligomers. Conclusions/Significance: Thermal denaturation of RIα leads to partial loss of native packing with exposure of aggregation-prone motifs, such as the B' helices in the phosphate-binding cassettes of both CNB domains. The topology of the β-sandwiches in these domains favors inter-molecular β-aggregation, which is suppressed in the ligand-bound states of RIα under physiological conditions. Moreover, our results reveal that the CNB domains persist as structural cores through heat-denaturation. © 2011 Dao et al

Ashwagandha Derived Withanone Targets TPX2-Aurora A Complex: Computational and Experimental Evidence to its Anticancer Activity

Cancer is largely marked by genetic instability. Specific inhibition of individual proteins or signalling pathways that regulate genetic stability during cell division thus hold a great potential for cancer therapy. The Aurora A kinase is a Ser/Thr kinase that plays a critical role during mitosis and cytokinesis and is found upregulated in several cancer types. It is functionally regulated by its interactions with TPX2, a candidate oncogene. Aurora A inhibitors have been proposed as anticancer drugs that work by blocking its ATP binding site. This site is common to other kinases and hence these inhibitors lack specificity for Aurora A inhibition in particular, thus advocating the need of some alternative inhibition route. Previously, we identified TPX2 as a cellular target for withanone that selectively kill cancer cells. By computational approach, we found here that withanone binds to TPX2-Aurora A complex. In experiment, withanone treatment to cancer cells indeed resulted in dissociation of TPX2-Aurora A complex and disruption of mitotic spindle apparatus proposing this as a mechanism of the anticancer activity of withanone. From docking analysis, non-formation/disruption of the active TPX2-Aurora A association complex could be discerned. Our MD simulation results suggesting the thermodynamic and structural stability of TPX2-Aurora A in complex with withanone further substantiates the binding. We report a computational rationale of the ability of naturally occurring withanone to alter the kinase signalling pathway in an ATP-independent manner and experimental evidence in which withanone cause inactivation of the TPX2-Aurora A complex. The study demonstrated that TPX2-Aurora A complex is a target of withanone, a potential natural anticancer drug

Preference of Small Molecules for Local Minimum Conformations when Binding to Proteins

It is well known that small molecules (ligands) do not necessarily adopt their lowest potential energy conformations when binding to proteins. Analyses of protein-bound ligand crystal structures have reportedly shown that many of them do not even adopt the conformations at local minima of their potential energy surfaces (local minimum conformations). The results of these analyses raise a concern regarding the validity of virtual screening methods that use ligands in local minimum conformations. Here we report a normal-mode-analysis (NMA) study of 100 crystal structures of protein-bound ligands. Our data show that the energy minimization of a ligand alone does not automatically stop at a local minimum conformation if the minimum of the potential energy surface is shallow, thus leading to the folding of the ligand. Furthermore, our data show that all 100 ligand conformations in their protein-bound ligand crystal structures are nearly identical to their local minimum conformations obtained from NMA-monitored energy minimization, suggesting that ligands prefer to adopt local minimum conformations when binding to proteins. These results both support virtual screening methods that use ligands in local minimum conformations and caution about possible adverse effect of excessive energy minimization when generating a database of ligand conformations for virtual screening

Normal-Mode-Analysis–Monitored Energy Minimization Procedure for Generating Small–Molecule Bound Conformations

The energy minimization of a small molecule alone does not automatically stop at a local minimum of the potential energy surface of the molecule if the minimum is shallow, thus leading to folding of the molecule and consequently hampering the generation of the bound conformation of a guest in the absence of its host. This questions the practicality of virtual screening methods that use conformations at local minima of their potential energy surfaces (local minimum conformations) as potential bound conformations. Here we report a normal-mode-analysis–monitored energy minimization (NEM) procedure that generates local minimum conformations as potential bound conformations. Of 22 selected guest–host complex crystal structures with guest structures possessing up to four rotatable bonds, all complexes were reproduced, with guest mass–weighted root mean square deviations of <1.0 Å, through docking with the NEM–generated guest local minimum conformations. An analysis of the potential energies of these local minimum conformations showed that 22 (100%), 18 (82%), 16 (73%), and 12 (55%) of the 22 guest bound conformations in the crystal structures had conformational strain energies of less than or equal to 3.8, 2.0, 0.6, and 0.0 kcal/mol, respectively. These results suggest that (1) the NEM procedure can generate small–molecule bound conformations, and (2) guests adopt low-strain–energy conformations for complexation, thus supporting the virtual screening methods that use local minimum conformations

Identification of Key Residues for pH Dependent Activation of Violaxanthin De-Epoxidase from Arabidopsis thaliana

Plants are often exposed to saturating light conditions, which can lead to oxidative stress. The carotenoid zeaxanthin, synthesized from violaxanthin by Violaxanthin De-Epoxidase (VDE) plays a major role in the protection from excess illumination. VDE activation is triggered by a pH reduction in the thylakoids lumen occurring under saturating light. In this work the mechanism of the VDE activation was investigated on a molecular level using multi conformer continuum electrostatic calculations, site directed mutagenesis and molecular dynamics. The pKa values of residues of the inactive VDE were determined to identify target residues that could be implicated in the activation. Five such target residues were investigated closer by site directed mutagenesis, whereas variants in four residues (D98, D117, H168 and D206) caused a reduction in enzymatic activity indicating a role in the activation of VDE while D86 mutants did not show any alteration. The analysis of the VDE sequence showed that the four putative activation residues are all conserved in plants but not in diatoms, explaining why VDE in these algae is already activated at higher pH. Molecular dynamics showed that the VDE structure was coherent at pH 7 with a low amount of water penetrating the hydrophobic barrel. Simulations carried out with the candidate residues locked into their protonated state showed instead an increased amount of water penetrating the barrel and the rupture of the H121–Y214 hydrogen bond at the end of the barrel, which is essential for VDE activation. These results suggest that VDE activation relies on a robust and redundant network, in which the four residues identified in this study play a major role

2- and 8-alkynyl-9-ethyladenines: Synthesis and biological activity at human and rat adenosine receptors

The synthesis of a series of 9-ethyladenine derivatives bearing alkynyl chains in 2- or 8-position was undertaken, based on the observation that replacement of the sugar moiety in adenosine derivatives with alkyl groups led to adenosine receptor antagonists. All the synthesized compounds were tested for their affinity at human and rat A1, A2A, and A3 adenosine receptors in binding assays; the activity at the human A2B receptor was determined in adenylyl cyclase experiments. Biological data showed that the 2-alkynyl derivatives possess good affinity and are slightly selective for the human A2A receptor. The same compounds tested on the rat A1 and A2A subtypes showed in general lower affinity for both receptors. On the other hand, the affinity of the 8-alkynyl derivatives at the human A1, A2A, and A2B receptors proved to be lower than that of the corresponding 2-alkynyl derivatives. On the contrary, the affinity of the same compounds for the human A3 receptor was improved, resulting in A3 selectivity. As in the case of the 2-alkynyl-substituted compounds, the 8-alkynyl derivatives showed decreased affinity for rat receptors. However, it is worthwhile to note that the 8-phenylethynyl-9-ethyladenine was the most active compound of the two series (Ki in the nanomolar range) at both the human and rat A3 subtype. Docking experiments of the 2- and 8-phenylethynyl-9-ethyladenines, at a rhodopsin-based homology model, gave a rational explanation of the preference of the human A3 receptor for the 8-substituted compound